Statistical Mechanics of an Urban Cross : A Solvable Model 



O 

o 
o 



(N 



I 

o 
o 



> 

\o 
cn 
i> 
o 
o 
o 



X3 
O 

o 



X 
J3 



M.Ebrahim Fouladvand^''^ , Masoud NcmatoUahi^ and Farzad Zolfaghari^ 
1 Department of Physics, Sharif University of Technology, P.O. Box 11365-9161, Tehran, Iran. 
2 Tehran Traffic Control Company (TTCC), P.O. Box 15836 , Tehran, Iran. 
3 Institute for Studies in Theoretical Physics and Mathematics ,P.O. Box 19395-5531, Tehran, Iran. 

(February 6, 2008) 

We propose a model for the intersection of two urban streets. The traffic status of the crossroads 
is controlled by a set of traffic lights which periodically switch to red and green with a total period 
ofT. Two different types of crossroads are discussed. The first one describes the intersection of two 
one-way streets, while the second type models the intersection of a two-way street with an one-way 
street. We assume that the vehicles approach the crossroads with constant rates in time which are 
taken as the model parameters. We optimize the traffic flow at the crossroads by minimizing the 
total waiting time of the vehicles per cycle of the traffic light. This leads to the determination of 
the optimum green-time allocated to each phase. 



Over the past decade, the majority of the studies in the 
inter-disciplinary field " physics of traffic " have been al- 
located to the highway traffic ( for a review see Ref. 0-^] 
) and less attention has been paid to the urban traffic. 
The simulation of traffic flow in a large-sized city is a 
formidable task and many degrees of freedom have to 
be involved (see e.g. [Q). In practice, streets of a city 
form a network of junctions that are linked togetether. 
Each crossroads receives demands (vehicles attempting 
to pass the cross ) and at each crossroads there exists a 
traffic light which, with some certain programming, con- 
trols the transportation. The first model introduced by 
statistical physicists for the description of the city net- 
work which is known as the BML model in the literature, 
uses a deterministic cellular automata framework j^] and 
predicts a sort of phase transition from free-flow to a 
jammed state, but the model assumptions are too sim- 
ple to be taken seriously for practical purposes. Later, 
a more serious model was introduced by Schadschneider 
and Chowdhury . This model combines the BML mode 
together with the Nagel-Schreckenberg model of highway 
traffic 0. Despite its nice formulation, there are still 
a lot of simplifications which prevent it from being an 
effective and applicable model to a city network. The 
empirical mechanism by which the traffic lights are con- 
trolled is generally divided into two distinct methods : 
fixed time and real-time. In the fixed-time method, a 
fixed value of time is allocated to the traffic light as well 
as its sub-phase times. In the real-time method, which is 
becoming increasingly popular in great cities, the ensem- 
ble of crossroads are intelligently controlled by a central 
controller. The control mechanism is usually based on 
the concept of producing a kind of green waves between 
the crossroads. These waves interact with each other, 
and, if the passing demands are of high values, the green 
waves may have destructive effects on each other, and, 
hence, the concept of green wave may fail to be the final 
solution for optimizing the overall flow. 
In this paper we aim to analyze a single crossroads in de- 



tail in order to find a better insight to the problem of op- 
timizing the total fiow in cities. Our single crossroads is 
the result of the intersection of two perpendicular streets. 
In their simplest structure, these streets can each direct 
a one-way traffic flow. With no loss of generality, we take 
them as one-way South to North (S-N) and West to East 
(W-E) streets. Cars arrive at the south and the west 
entrances of the crossroads. Here we assume the arrival 
rates of the cars , i.e. the number of cars reaching the 
crossroads per second, are constant in time. Although 
from everyday driving experiences in cities, we know that 
these rates have fluctuations in the course of time; yet in 
definite time intervals , the assumption of the constant 
arrival rates could be justified at least on a average level. 
We take the arrival rates to be ai (for S-N cars ) and 
a2 (for the W-E cars) respectively. Also we denote the 
passing-rate of cars ( number of cars passing the cross- 
roads in the unit of time during the green-phase ) by (3i 
and /?2- The period of the traffic lights is taken to be a 
definite value T which is assumed to remain constant in 
the time interval under consideration. The starting time 
of each cycle of the traffic light is at the moment which 
the light turns green for the S-N street. The S-N light 
remains green until Ti. At Ti the traffic lights turn red 
for the S-N street and simultaneously changes to green 
for the W-E street. This is the beginning of the second 
phase which continues from Ti to T (end of one cycle). 
During Phase I ( < i < Ti ) , the S-N cars can pass the 
crossroads northward and W-E cars are stopped for the 
red fight. Over Phase II ( Ti < i < T ), the S-N cars 
must wait behind the red light whilst the W-E cars are 
passing through the crossroads. 

Now the basic question is how should traffic engineers 
adjust the value of Ti in order to optimize the traffic fiow 
throughout the intersection? Regarding the fact that 
drivers are identical human beings and with the assump- 
tion that the number of passengers of each car takes an 
equal average value for each direction, the optimization 
task is realized by minimizing the total waiting time of 
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the cars per cycle of the traffic lights. For this purpose, 
we introduce two quantities A^i and which represent 
the number of cars stopping at the red lights in the red 
phases of the S-N and W-E streets respectively. Clearly 
A^i and A^2 are functions of time and, in general, are di- 
vided into different lanes on the streets. For the first 
cycle of the steady time-interval ( in which the arrival 
rate are taken to be constant ) the following equations 
are easily written down: 



Ari(Ti) = [7Vi(0) + (ai-/3i)Ti]^ 



7V2(Ti) = 7V2(0)e + a2Ti 



(1) 



(2) 



iVi(T) = [iVi(O) + (ai - /3i)ri]0 + ai(r - T^) (3) 



N2{T) = [N2{0)e + a2T - /32(T - Ti)]e 



(4) 



The symbols ensures the positiveness of the quan- 
tities in the brackets, i.e. the value of is one if the 
quantity in the bracket is positive and zero elsewhere. 
This limitation is dictated to us since, by definition, the 
quantities Ni and N2 can only take positive values. 

We now define the Total Waiting Time (TWT) of the 
vehicles per cycle of the traffic lights. It is the total time 
wasted by the vehicles during their stop in the red phases. 
The TWT is sum of the sub-waiting times of each direc- 
tion. Denoting the TWT and the sub waiting times by 
T^, r(^'i) and r(^'2) respectively, the following equa- 
tions could be written for the n-th cycle 



rT.{W,l) ^ ^ irp 



Tif + Ni{nT + Ti){T-Ti) (5) 



(6) 



In the above equations, one must take into account the 
fact that during the red phase of a street, not all the cars 
reaching the crossroads wait equally. The cars reach- 
ing the intersection just after the light goes red ( un- 
lucky drivers) have the most waiting time, while those 
reaching the crossroads just before the lights turns green 
(lucky drivers) spend no time waiting. This leads to an 
average factor ^ in the above equations. The value of 
N2{Q) refers to the number of the stopped cars in the W- 
E street just after the instant when the light turns red for 
the W-E cars (in the first cycle). Here two distinct situa- 
tions are identified. In the first case, N2 (0) is nearly zero 
which describes the situation when the W-E street is not 
crowded. The other case is N2{Q) > which corresponds 
to the situation of heavy traffic in the W-E street. This 
means that even in the final stages of the green phase for 
the W-E direction, still there is a considerable number of 
W-E cars desiring to pass the cross eastwards. The same 
scenario goes with the value of A''i(Ti), and, as we will see 



in what follows, these different cases lead to different pre- 
scriptions of the optimization. Before investigating these 
cases in details, it should be mentioned that for the ar- 
bitrary values of the rates we do not have a periodical 
structure in the values of the Nx and A^2- This can be 
easily verified by comparing Ni{T) and N2{T) with their 
values at i = which implies that, in general, one has 
Ni{0) ^ N^{T) for z = 1,2. Now let us discuss the situ- 
ation where A^i(Ti) = N2{0) = 0. This refers to a light 
traffic flow in both directions during the first cycle of the 
interval. One can easily check that if the following stabil- 
ity conditions aaT - /?2(T - Ti) < and aiT - PiTi < 
hold. Then we always have the similar conditions in the 
next cycles Ni{nT + Ti) = N2{nT) = where n denotes 
the number of the traffic light cycle. In this particular 
case which here after will be called State I, the TWT 
is independent of the cycle-nmnber. and minimizing the 
TWT with respect to Ti leads to the following equation: 



ai 



ai + a2 



(7) 



this answer must be consistent with the stability condi- 
tions which leads to the constraints /32 > cti + ct2 and 
Pi > ai+ a2 among the rates. We now investigate a to- 
tally different situation i.e. a crowded crossroads. Sup- 
posing that the first cycle is characterized by the con- 
ditions N2{0),Ni{Ti) > 0. one could easily verify that 
provided the relations ai > (3i and a2T > (52{T — T\) 
hold, we have a quasi-stationary condition in the next 
cycles i.e. N-i{nT + Ti), N2{nT) > 0. In sharp contrast 
to the State I, (this state which will be referred to as the 
State II), the values of iVi and N2 would be functions of 
the cycle number. It can be shown that in the large cycle- 
number limit one has T^^^\_^ ~ n(T - ri)(aiT - PiTi) 

and r^^'^^j^ ^ nTi{a2T — f32{T — Ti)) . We now minimize 
the TWT with respect to Ti which leads to the following 
equation: 



Ti = 



f3i+ 02 + ai-C(2 r 
2(/3i+/?2) 



(8) 



the consistency of this solution with stability conditions 
yields the following constraints: ai > /3i and ail32 + 
2a2(3i + a2p2 > /9i/32 + Pi- -^Iso positiveness of Ti itself 
imposes the extra restriction ai — 02 < Pi +02- Next we 
consider the situation (State III) where in the first cycle 
of the traffic fight, one has Ni{Ti) = but ^¥2(0) > this 
correspond to the situation in which the S-N street has 
a light traffic flow while the W-E street has a heavy one. 
The conditions for having these situations repeated in the 
next cycles are aiT-PiTi < and a2T-P2{T-Ti) > 0. 
In the large cycle- number limit, the minimizing of the 
TWT leads to the value 



Ti = 



P2 



a2 



2p2 



T 



(9) 



It could be easily checked that the above solution is in- 
consistent with the stability conditions, and, hence, is not 
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acceptable as an optimization of traffic light. This incon- 
sistency will remain in the final state ( State IV ) which 
is characterized by Ni{nT + Ti) > and 7V2(nT) = 0. In 
this case, the solution obtained by minimizing the TWT 
is again inconsistent with the stability conditions. In 
these last two states, one must numerically find the con- 
sistent minima. We next consider another more usual 
configuration of a crossroads. Here we let the vehicles 
move in the both directions on the S-N street but still 
the vehicles in the W-E street are restricted to move 
eastwards. This situation describes a one-way to two- 
way urban intersection. Here each cycle consists of three 
phases. In the first phase which lasts for < i < Ti, the 
traffic light is green for the S-N cars and red for the other 
two directions. During the second phase which starts at 
Ti and finishes at T2, the traffic light is green for the N- 
S cars and rod for the other two directions. In the final 
phase, which last for T2 <t <T, the traffic light remains 
green for the W-E cars and red for the N-S as well as 8- 
N directions. The entrance rates are taken to be ai,a2 
and as and we denote the passing rate by Pi , P2 and P3 
for each direction respectively. The starting time of the 
cycles is chosen to be the moment at which the traffic 
light turns green for the S-N direction. Similar equations 
for the Vehicle-Numbers Ni , N2 and A^3 could be written 
down and one in principle can evaluate the TWT in terms 
of these quantities. The exact form of the TWT strongly 
depends on the positiveness of the vehicle-numbers just 
after the traffic lights goes red for the respective direc- 
tion. In the case under consideration, eight different pos- 
sibilities can occur. Let us only discuss the most proba- 
ble one which corresponds to the the following situation: 
Ni{nT + Ti) = N2{nT + T2) = Ns{nT) = In this case 
the traffic status is light in each of the three directions. 
Following the same discussion made for the previous 
crossroads, it could be easily verified that the triple sta- 
bility condition for the validity of the above assumptions 
arc as follows: aiT - PiTi < 0, a2T - (32iT2 - Ti) < 
and Q'sT — (3s{T — T2). It can be shown that minimizing 
the TWT with respect to Ti and T2 
leads to the following fixation of Ti, T2 



Ti=T 



To = T- 



a\{a.2 + as) — 0:20:3 
q:i(q;2 + OLz) + 0:203 

20102 
01(0:2 + 03) -I- 0203 



(10) 



(11) 



One directly observes that in the symmetric case of 
equal arrival rates, T\ and T2 take the expected values 
^ and ^ respectively. The other point which must be 
mentioned is that in this traffic status, T\ and T2 do not 
depend on the passing rates, and are solely determined 
by the arrival rates. The remaining seven status will be 
discussed elsewhere. 

For comparison of our model to the empirical data, a 
time-series analysis on two of Tehran intersections were 
carried out. The central part of Tehran is under control 



of the SCATS [Sydney Coordinated Adaptive Traffic Sys- 
tem) that is an intelligent traffic controller system. The 
strategy followed by most of urban traffic control sys- 
tems is based on establishing green-wave along the ma- 
jor streets of cities. We considered two diflPerent intersec- 
tions. The first one which is located in Tehran downtown 
connects Valiasr St. to Takhtejamshid St. Both of these 
are one-way and major streets. The other crossroads con- 
nects Abbasabad St. (one-way) to Mahnaz St. (one-way) 
.In contrast to the previous case, here the first street is a 
major while the second street is a minor one. The data 
set is provided by magnetic counting loops which are in- 
stalled just before the pedestrian-line of each crossroads. 
The data were collected on second of July, 2000. The 
three following figures show a time-series of data for our 
major-to-minor crossroads. 




Fig. 1: The cycle-ratio of green-time of Abbasabad street 
to the total period of the traffic light cycle. The numbers on 
the horizontal axes are the cycle number. 



Fig. 2: The cycle-ratio of the passed vehicles (during the 
green phases) of Abassabad street to Mahnaz street of each 
traffic light cycle. The numbers on the horizontal axis denote 
the cycle number. 
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Fig. 3: the above graph denotes the total number of passed 
vehicles (from both streets) during fifteen-minute time inter- 
vals. The lower graph shows the number of passed vehicles 
form Mahnaz street during fifteen-minute intervals. 
As seen from the graphs, the cycle-ratio of the green time 
allocated to Abbasabad str. strongly fluctuates due to 
demand fluctuations received by the cross. In order to 
have a rough estimation of our model parameters, we 
considered a two-hour time interval between 12:30-14:30 
. We approximated the ratio of ^ by the ration of the 
total number of vehicles which have passed through Ab- 
basabad Str. to that number for Mahnaz Street. This 
yields the value ^ = 0.35. Putting this value into equa- 
tion (7) yields ^ = 0.74. On the other hand the averaged 
value of the cycle-ratio of the green time of Abbasabad 
street over the two-hour period leads to result ^ = 0.64. 
The following flgures belong to our major-to-major cross- 
roads. 



i 09 




s 

J 0.7 
1 

1 OB 




1 05 




£ 




% 0.1 




100 200 300 

7:30 10:30 13:30 16:30 19:30 22:30 



Fig. 4: the cycle-ratio of green time allocated to Takhte- 
jamshid Street. The numbers on the horizontal axis denote 
the cycle number. 
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Fig. 5: cycle-ratio of the passed vehicles from Takhte- 
jamshid street to Valiasr street. 
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Fig. 6: the above graph denotes the total number of passed 
vehicles (from both streets) during fifteen-minute time inter- 
vals. The lower graph shows the number of passed vehicles 
form Takhtejamshid str. during fifteen-minute intervals. 
Here we focused on the interval 13:30 - 15:30. The aver- 
aged value of ^ is 0.52 while the same amount evaluated 
from eq.(7) is 0.56. 

In conclusion, we have developed a prescription for the 
traffic-light programming at a single urban crossroads. 
The method is based on minimizing the TWT of cars 
stopping in the red phases of the traffic light. The value 
of T should be so tuned that during the green phases, 
time-headway is less than a certain value. In fact, in the 
model the values of the passing rate refer to the maxi- 
mum capacity of cross ( maximum number of cars passing 
the cross in the unit of time ), which is plausible if the 
T is appropriately adjusted with the congestion of the 
crossroads. Optimizing the traffic at each crossroads is 
the stating point of the more comprehensive problem of 
city network. Nevertheless, our model best suits those 
marginal intersections of cities where the effect of the 
other crossroads is suppressed. 
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